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Abstract 



The NA48 experiment at CERN has performed a new measurement of direct CP 
violation, based on data taken in 1997 by simuhaneously cohecting Kl and Kg 
decays into and ■k'^tt . The resuh for the CP violating parameter Re (e'/e) is 
(18.5 ± 4.5(stat) ±5.8(syst)) x 10"^. 
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1 Introduction 

CP violation has been observed so far only in the decays of neutral kaons [|l| . The 
main component of the effect |0] occurs in the mixing between and eigenstates. 
The physical states Ks and Kl deviate from pure CP=±1 eigenstates, with the mixing 
described by the parameter e. On the other hand, direct CP violation can occur in the 
transition from the neutral kaon eigenstate with CP=— 1 to a two-pion final state, which 
has CP=+1. This can be observed through the interference of the amplitudes of different 
final state isospin and is described by the parameter e'. The measured quantity is the 
double ratio of decay widths: 

The Standard Model relates CP violation to the existence of three generations of 
quarks and to the complex phase present in the CKM matrix [Q. Computations based 
on the Standard Model typically predict Iie{e'/e) in the range to 10 x 10~^ The 
first evidence for direct CP violation was published in 1988 Subsequent measurements 
(23.0 ± 6.5) X 10~^ and (7.4 ± 5.9) x lO^'' [§ have been only marginally consistent 
with each other; the measurement in |^ shows evidence for direct CP violation while 
that in is consistent with no effect. The KTeV collaboration has recently measured 
(28.0 ±4.1) X 10~^ using a technique similar to the one used in |p. 

In this letter, a new measurement of Re [e'/e) obtained by the NA48 collaboration 
at CERN is presented, based on the first data collected with the experiment in 1997. 
While the ultimate accuracy of the NA48 experiment will be significantly better than 
that reported in this study, this result is of sufficient precision to help clarify the current 
experimental situation. 

2 Principle of measurement 

The measurement of Re {e'/ e) proceeds by counting the number of events in each 
of the four decay modes and computing the double ratio R. The NA48 experiment is 
designed in such a way that many systematic uncertainties can be controlled with high 
accuracy by exploiting cancellations of effects which contribute symmetrically to different 
components of the double ratio. 

Data are collected simultaneously in the four decay channels, minimizing the sensi- 
tivity of the measurement to accidental activity and to variations in beam intensity and 
detection efficiency. Kl and Ks decays are provided by two nearly collinear beams, with 
similar momentum spectra, converging in the centre of the main detector. 

In order to reduce the difference in acceptance due to the large difference in average 
decay lengths, only K^ decays occurring in the region also populated by Ks decays are 
used for the measurement of the double ratio. Furthermore, a cancellation of the residual 
difference is obtained by weighting the Kl events used in the double ratio with a function 
of the proper lifetime r, proportional to the expected ratio of Ks and Kl decays at time 
r. In this way, the systematic accuracy of the result does not rely on a detailed Monte 
Carlo simulation of the experiment. 

3 Beam and detector 
3.1 The kaon beams 



The neutral beams ||T^ are derived from 450 GeV/c protons extracted from the 



CERN SPS. Due to the different mean decay lengths (Al=3.4 km, As'=5.9 m at the average 
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kaon momentum of 110 GeV/c), two separate production targets are used, located 126 m 
and 6 m upstream of the beginning of the decay region. 

For each SPS pulse (2.4 s spill every 14.4 s), 1.1 x 10^^ protons hit the pro- 
duction target. A neutral secondary beam, with ±0.15 mrad divergence is generated at a 
production angle of 2.4 mrad. Three stages of collimation are used with the final element 
located just upstream of the beginning of the decay region. 

Part of the non-interacting primary protons impinge on a bent silicon mono-crystal 



Tl| . A small fraction of these protons undergo channelling, and are bent to produce a 
coUimated beam of 3 x 10^ protons per pulse transported to the Ks production target. 

The Ks collimator selects secondary neutral particles at 4.2 mrad production angle, 
with a divergence of ±0.375 mrad. The different production angles for the Kl and Ks 
beams are chosen to minimize the difference in decay momentum spectra. The Ks beam 
enters the fiducial decay volume 68 mm above the Kl beam. The beams converge with an 
angle of 0.6 mrad and the axes of the two beams cross at the position of the electromagnetic 
calorimeter. The total flux of Kl through the fiducial region is ~ 1.5 x lO'' per spill, and 
the Ks flux is ^ 2 X 10^ per spill. 

The protons directed to the Ks target are detected by an array of scintillation coun- 



ters which comprise the tagging detector [0. This device is used to tag Ks decays. The 



time of any protons and the time of the event in the main detector are each reconstructed 
relative to a common free running clock. The presence (absence) of a proton reconstructed 
in coincidence with the event defines the event as Ks {Kl). The signals are digitized by 
960 MHz FADCs. The proton time resolution is 120 ps and the double-pulse separation 
is 4 ns. 

The neutral beam produced at the Ks target traverses an anti-counter (AKS), 
formed by a set of three scintillation counters preceded by an aligned 3 mm thick iridium 
crystal which enhances photon conversions (1.8 effective radiation lengths) [jl3|. This de- 
tector is used to veto Ks decays occurring upstream, thus providing an accurate definition 
of the beginning of the decay region. 

The decay region is contained in an evacuated (< 3 x 10~^ mbar) 90 m long tank 
with a 0.9 mm (0.003 radiation lengths) thick polyimide (Kevlar) composite window at 
the end. The window extends from an inner diameter of 18.4 cm to an outer diameter of 
2.3 m. The neutral beam continues in a 16 cm diameter evacuated tube to a {Kl) beam 
monitor and the beam dump, both downstream of all the detector elements. 

3.2 Main detector 

The layout of the main detector is shown in figure |[ Kaons which decay into charged 
pions are measured by a magnetic spectrometer comprised of four drift chambers [|T1|, |TB|, 



and a dipole magnet giving a horizontal transverse momentum kick of 265 MeV/c. The 
volume between the drift chambers is filled with helium at atmospheric pressure. The 
drift chambers are each comprised of eight planes of sense wires, two horizontal, two 
vertical and two along each of the 45° directions. Only horizontal and vertical planes 
are instrumented in DCH3. The resolution in each projection is 90 /im and the average 
efficiency is 99.5 % per plane. The geometrical accuracy is better than 0.1 mm/m. The 
momentum resolution is (Tp/p ~ 0.5 % © 0.009 % x p (p in GeV/c), where © means that 
the contributions should be added in quadrature. 

Two plastic scintillator hodoscope planes are placed after the helium tank. They 
are segmented in horizontal and vertical strips, respectively, with widths varying between 
6.5 cm and 10 cm, arranged in 4 quadrants. Fast logic provides a signal for the first level 
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of the 7r"'"7r~ trigger. This trigger requires the presence of tracks in two opposite quadrants 
of the hodoscope. The hodoscope data are also used in the analysis to determine the event 
time of TT+TT" events for the tagging procedure. The time resolution is ~ 200 ps per track. 

A quasi-homogeneous liquid krypton electromagnetic (e.m.) calorimeter with 127 cm 
(27 radiation lengths) long projective tower readout is used to measure the photons from 
7r°7r° events. The electrode structure consists of copper/beryllium ribbons extending be- 
tween the front and the back of the detector with a ± 48 mrad accordion geometry. The 
13212 readout cells each have a cross section of 2 x 2 cm^ at the back of the active region. 
The projective towers point to the average Ks decay position (110 m upstream of the 
detector) so that the measurement of the photon positions is insensitive to the initial 
conversion depth. The transverse scale at the back plane of the calorimeter is measured 
to better than 0.2 mm/m. The initial current induced on the electrodes by the drift of the 
ionization is measured using pulse shapers with 80 ns FWHM and digitized with 40 MHz 
FADCs The energy resolution is a{E)/E ~ 0.125/E©0.032//E©0.005 {E in GeV). 
The spatial resolution is better than 1.3 mm and the time resolution is better than 300 ps 
for 20 GeV photons. 

The e.m. calorimeter incorporates a detector made of a 4 mm thick plane of scintil- 
lating fibres located near the maximum of the shower development. This device provides 
an independent trigger for vr^vr'' decays, and an additional measurement of the event time. 

An iron-scintillator hadron calorimeter, 6.7 nuclear interaction lengths thick, is used 
to complement the krypton calorimeter to measure the energies of hadrons in the trigger. 
Three planes of scintillation counters, segmented into 25 cm wide counters, shielded by 
80 cm thick iron walls, are used for muon identification to reduce the background from 
Kl — > TT/iz/ (A'^s) events. Seven counter arrays (indicated as ANTI 1 - ANTI 7 on figure p, 
consisting of plastic scintillators and iron converters, surround the vacuum and helium 
tanks to detect photons which miss the calorimeter. In the present analysis, they are only 
used to monitor accidental activity. 
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Figure 1: Layout of the main detector components. 
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3.3 Trigger and DAQ 

The readout of the e.m. calorimeter sends the calibrated analogue sums of 2 x 8 cells 
to the trigger. The trigger digitizes the data and operates as a 40 MHz 



pipelined computation system with a latency of 3 fis. The calorimeter data are reduced 
to X and y projections, each consisting of 64 strips, which are used to reconstruct the 
total deposited energy £"0, the first and second moments of the energy distribution and to 
reconstruct and measure the time of energy peaks. This information is used to reconstruct 
the radial position of the centre of gravity C of the event and the proper decay time t/ts 
in units of the mean Ks lifetime ts- Events are selected if Eo > 50 GeV, C < 15 cm, 
t/ts < 5.5 and if there are no more than five peaks in both views in a 9 ns time window. 
The ir^n^ trigger rate is about 4500 per pulse and the deadtime is negligible. 

A two level tt+tt" trigger is used. At the first level, the hodoscope trigger is placed 
in coincidence with a total energy condition (> 30 GeV), defined by adding the hadron 
calorimeter energy to Eq from the trigger. The first level trigger is downscaled by a 
factor of two, the resulting trigger rate is about 100,000 per pulse. The second level 7r~^7i~ 



trigger [|18l uses information from the drift chambers. It consists of hardware coordinate 
builders for each view and a farm of asynchronous microprocessors which reconstructs 
tracks. Triggers are selected if the tracks converge to within 5 cm, their opening angle 
is less than 15 mrad, the reconstructed proper decay time is less than 4.5 ts and the 
reconstructed kaon mass is at least 95 % of the nominal mass. The maximum decision 
latency is 100 /xs. The second level trigger rate is about 2000 per pulse. 

Downscaled triggers, based on relaxed conditions are collected to compute the effi- 
ciency of the main triggers. "Random events" record any activity in the detectors 69 //s 
after a heavily downscaled beam counter signal. The 69 /is delay corresponds to the pe- 
riodicity of the slow proton extraction (three revolution periods) of the SPS accelerator. 
These events contain signals which are representative of accidental activity in the detec- 
tors. 

The overall deadtime of the trigger system is below 1 %, the largest contribution 
being due to the tt+tt" trigger. The deadtime condition is recorded and made completely 
symmetric offiine. The detector readout systems provide information on all the activity 
in the detectors in a readout window of about 250 ns surrounding the trigger time. The 
readout window for the drift chambers is larger, about 1.2 /is. 

4 Event reconstruction and selection 
4.1 7r°7r° mode 

Photon showers are found in the calorimeter by looking for maxima in the digitized 
pulses from individual cells in both space and time and accumulating the energy within a 
radius of 11 cm. Photon energies are corrected for: energy outside the cluster boundary; 
energy lost in non- working cells (about 0.4% of the channels); energy from the other 
photons in the event; small variations of the energy measurement depending on the impact 
point within the cell and small space charge effects from the accumulation of positive 
ions [^] because the calorimeter was operated at half the nominal 3kV drift voltage in 
1997. 

K — > 7r°7r° decays are selected by requiring four clusters with energies between 
3 GeV and 100 GeV which are in time within ±5 ns of the average of the four. Fiducial 
cuts are applied to ensure that photon energies are well measured. These cuts include 
the removal of photons within 2 cm of a non- working cell and photons within a ± 4 cm 
wide column due to a failure in a HV connection during the 1997 run. The minimum 
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distance between photon candidates is required to be greater than 10 cm. To reduce the 
background from Kl ^ 3 tt" decays, events with an additional cluster of energy above 1.5 
GeV and within ± 3 ns of the tt'^tt" candidate are rejected. 

The kaon energy is obtained with a resolution of ~ 0.6 % from the sum of the 
photon energies. The longitudinal position of the decay vertex relative to the front of the 
calorimeter D is reconstructed from the energies Ei and positions Xj, Ui of the four selected 
clusters, assuming that their invariant mass is the kaon mass (m/<) as follows: 



V * j>« 



(2) 



and is used to reconstruct the proper decay time r. The resolution of D is 50 to 70 cm 
depending on energy. The invariant masses mi and m2 of the two photon pairs are com- 



puted using D and compared to the nominal vr" mass 
computed as follows: 



For this, 



a X 



variable is 



X 



(mi + m2)/2 - m^o 
0-+ 



+ 



(mi - m2)/2 



(3) 



where (t± are the resolutions of (mi ± m2)/2 observed in the data (parameterised as a 
function of the lowest photon energy). The three-fold choice of photon pairings is made by 
selecting the combination with the lowest x^. Due to the constraint on m^^, (mi ±m2) are 
to good approximation uncorrelated. Typical values of 0"+ and (T_ are 0.45 and 1.1 MeV/c^. 
To select good ti^ti^ candidates, the best combination is required to have < 13.5. This 
requirement removes about 7 % of Ks — > tt^tt'^ events in a pure Ks beam where there is no 
background. About 70 % of these are due to vr" Dalitz decays and to photon conversions in 
the Kevlar window or in the spectrometer volume upstream of the magnet. The remaining 
part is due to the e.m. calorimeter response, and is dominated by occasional photon- 
nucleon interactions, which result in a reduction of the observed energy. Both effects are 
highly symmetric between Ks and K^, and do not affect the measurement of the double 
ratio. 

The event time as used in the tagging procedure for tt'^tt^ events is computed com- 
bining eight time estimators from the two most energetic cells of each cluster. 



4.2 tt+tt" mode 

Tracks are reconstructed from the hits and drift times in the spectrometer. A vertex 
is defined at the point of closest approach of two oppositely charged tracks (provided 
they approach to within 3 cm). The momenta of the tracks are derived using a detailed 
magnetic field map. The energy of the event is reconstructed using the angle between the 
tracks before the magnet and the ratio of the track momenta, assuming the two tracks 
originate from K tt^ti" . This avoids systematic uncertainties from the knowledge of 
the magnetic field in the definition of the fiducial energy cut and lifetime measurement. 
Small corrections on the tracks are made for the field (integral ~ 2 x 10~^ T m) in the 
decay region before the spectrometer. 

The vertex resolution is typically 50 cm in the longitudinal direction z, and 2 mm in 
the transverse directions x, y. Since the beams are separated by about 6 cm in the decay 
region, a clean identification of Ks tt+tt^ and K^ tt+tt" decays is possible (this is 
referred to as "y- vertex tagging" in the following). 
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A variable A related to the decay orientation in the kaon rest frame is defined from 
the two track momenta pi and p2 as A = \pi — P2\/{pi + ^2)- Kaon energy dependent 
cuts A < 0.62 and A < 1.08 — 0.0052 x (where is the kaon energy in GeV) are 
made in order to avoid using events with tracks near the beam hole for which a precise 
Monte Carlo modelling is required. This cut also entirely removes A and A decays which 
are present in the Ks beam. 

To reject background from semileptonic Kl decays, tracks consistent with being 
either an electron or a muon are rejected. To identify electrons, the ratio E/p of the 
energy of a matching cluster in the e.m. calorimeter to the track momentum is computed. 
Both tracks must satisfy E/p < 0.8. This requirement reduces the Kl — > vrez/ (-ft'es) 
background by a factor of 500 while removing about 5 % of the 7r~^7i~ signal. Tracks are 
identified as muons if hits are found in time in the muon counters near the extrapolated 
track impact point. Events with identified muons are rejected. The rejection against i^^^a 
background is about 500 with an efficiency for 7r"^7r~ of 97 % (the loss coming essentially 
from TT^ decays to /i^). 

For good 7r~^TT~ events, the invariant mass m^^r should be equal to the kaon mass. The 
invariant mass resolution am is typically 2.5 MeV/c^. An energy dependent cut at ±3(Tm 
is applied. A further reduction of background due to semileptonic decays is obtained 
using a cut based on the kaon transverse momentum. To define a selection which is 
symmetric between Ks and K^ decays, a variable p/, defined as the component of the 
kaon momentum orthogonal to the line joining the production target (identified with y- 
vertex tagging) and the point where the kaon trajectory crosses the plane of the first drift 
chamber, is used. The p/ resolution is independent of the target from which the event 
originated. To select n'^n' candidates, the cut pt'"^ < 2 x 10""^ GeV^/c^ is applied. Event 
losses (~ 8 X 10~^) are mainly due to large deviations in the vertex reconstruction which 
are independent of whether the event is K^ or Ks- 

The event time as used in the tagging procedure in the tt+tt^ mode is measured by 
combining the four times from the two hodoscope planes for each track. The procedure 
which is applied protects against tails in single time measurements (predominantly due to 
accidental hits). Events for which this procedure cannot be applied are rejected, removing 
2.3 X 10^^ of the 7i~^7i signal. 

4.3 Tagging 

The event is tagged by classifying it as Ks if there is a proton within a ± 2 ns 
coincidence window centred on the event time, otherwise the event is classified as Kl as 
shown in figure @a,c. 

An accidental coincidence between a proton traversing the tagger and a Kl event 
may cause the event to be wrongly counted as Ks- This probability is measured for 
7r~^n~ events to be = 0.1119 ± 0.0003 using ?/- vertex tagging. The measurement of 
the double ratio is sensitive to the difference in this probability between 7r°7r° and vr+vr" 
events, which is measured to be af/s ~ oi~ls ~ (^0 ± 5) x 10~^ leading to a correction to 
R of (18 ± 9) X 10"'^. This is measured by looking at proton coincidences at times which 
are offset from the actual event. Differences in accidental rate between event time and 
offset coincidence intervals are checked with Kl — >■ tt+tt" events, identified from y-vertex 
tagging, and with Kl — »■ 37r*^ events. The observed difference in effective rate seen by 
7r°7r° and tt+tt" events is compatible with the intensity dependence of the ti^ti~ trigger 
efficiency. 

Inefficiencies in reconstructing a proton or errors in measuring the event time or the 
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proton time can cause a coincidence to be lost and a true Ks event would be classified as 
Kl. From the vertex tagged ti^ti~ events, the probability for this to happen is a~gY, = 
(1.5±0.1) X 10~^ (figure^). Detailed studies show that this is dominated by inefficiencies 
in the tagger which are intrinsically symmetric between Tr^Tr*^ and tt+tt" decays. A sample 
of neutral events containing vr*^ Dalitz decays and photon conversions is used to compare 
the time measured using the calorimeter with the time of the tracks measured using the 
hodoscope (figure |^). This comparison provides the difference between tagging efficiencies 
for vr^vr'^ and tx^tx^ events, which is found to be zero, with an accuracy of ± 1 x 10~^, 
leading to an uncertainty on i? of ±6 x 10~^. As an independent check, Tr'^vr^ events 
containing a Dalitz decay have been used to y-vertex tag TX\^ events leading to a direct 
measurement of a^g^^ = (2.3l5;o) ^ 10~^, in agreement with the method above. 
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Figure 2: (a) Distribution of the minimum difference between tagger and event times for 
7r+7r~ decays. The peak corresponds to Ks events where there is a coincidence, (b) The 
same distribution, separated into Ks and K^ events using y-vertex tagging. This shows 
that the tails in the coincidence are very small, (c) as (a) for vr^vr'^ mode, (d) Coincidence 
time for neutral events with tracks (see text). 



4.4 Data quality selections, data samples 

Data collected in the vr+vr" mode are affected by an overflow condition in the drift 
chambers which resets the front-end readout buffers when there are more than seven hits 
in a plane within a 100 ns time interval [0. This occurs when an accidental particle 
generates an electromagnetic shower upstream of the spectrometer and sprays the drift 
chambers with particles. Overflows in each plane are recorded for at least 400 ns before 
and after each event. To maintain the highest detector efficiency, we require that events 
used in this analysis do not have overflows within ±312 ns of the trigger time. This 
interval, which is larger than the sum of the maximum drift time and reset time, excludes 
any correlation between overflows and selected events. To avoid indirect biases related to 
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different average beam intensities, the selection on overffows is also applied to the 7r°7r° 
mode resulting in a loss of 20 % of the events. 

The fiducial range in kaon energy and proper time used to count events is 70 < 
Ek < 170 GeV and < r < 3.5 ts, where r = is defined at the position of the AKS 
counter. For Kl, the decay time cut is applied on reconstructed r, while for Ks the cut 
at r = is applied using the AKS to veto events occurring upstream. The nominal r = 
position differs by 21.0 ± 0.5 mm between vr'^vr'' and tt'^tt" modes, which is accounted for 
in the analysis. The veto efficiency of the AKS is 0.9964 ± 0.0003 (0.9923 ± 0.0010) for 
TT^ir^ (vr+vr") decays and leads to a negligible correction to R. 

All events are required to satisfy a symmetric cut on the position of the centre of 
gravity. For vr^vr^ decays, the centre of gravity is defined as the energy-weighted average 
X, y position of the four showers at the face of the e.m. calorimeter. For tt+tt" decays it is 
the momentum-weighted average position of the track segments, measured upstream of 
the spectrometer magnet and projected onto the face of the e.m. calorimeter. The radial 
position of the centre of gravity is required to be within 10 cm of the beam axes. The 
sensitivity to differences in beam halo and to Ks scattering in the collimator and in the 
AKS counter is minimized by applying the same cut to both modes. For scattered Ks 
events, this cut is stronger than the p/ cut applied in the selection of charged events. 
It is however relatively wide compared to both the K^ and Ks beam-spot radii, which 
are 3.6 cm and 4.6 cm respectively, so that effects related to resolution smearing can be 
neglected. 

Table |l] shows the number of events collected in each channel after removing back- 
ground and correcting for mistagged events. The numbers are shown here without lifetime 
weighting. 



Table 1: Statistical samples (in thousands of events). 



Kl ^ n'n' 


Ks ^ vt^tt'^ 


Kl TT^TT 


Ks vr"^7r 


489 


975 


1,071 


2,087 



5 Computation of R and systematic uncertainties 

In order to be insensitive to the difference between the beam spectra (about ± 10 % 
in the range 70-170 GeV), the computation of R is done in 20 energy bins. The numbers of 
Ks and weighted Kl candidates are corrected for the tagging accidentals and inefficiencies 
as described above. Corrections for trigger efficiency, background subtraction, and residual 
acceptance difference between Ks and Kl are applied separately in each energy bin. The 
average of R over the bins is done with an unbiased estimator of log(i?). 

5.1 Trigger efficiency 

The efficiency of the vr'^vr" trigger is determined from events triggered by the scin- 
tillating fibre detector in the e.m. calorimeter. The average vr^vr^ trigger efficiency is mea- 
sured to be 0.9988 ± 0.0004 and is the same for Ks and Kl- 

The efficiency of the 7r"*"7r~ trigger is measured on a selected sample of events from an 
auxiliary trigger which reproduces the first level logic. The combined efficiency of the first 
level coincidence and the second level is 0.9168 ib 0.0009. The major contribution to the 
inefficiency comes from a timing misalignment in the trigger logic which existed for part 
of the run. Additional inefficiencies in the hodoscope trigger amount to (68 ± 8) x 10"^. 
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The double ratio is corrected by (+ 9 ± 23) x 10 ^ for the difference between the tt^tt 
trigger efficiency for Ks and r-weighted Kl decays. 

5.2 7r°7r° background 

The background in the Kl tt^tt^ sample originates from Kl Stt'' decays with 
two undetected photons. 

The estimation of the residual background after all cuts is done by extrapolating the 
number of events in the control region 36-135 to the signal region < 13.5 as shown 
in figure 0. The extrapolation factor for the background is measured using a sample of 
Kl Stt^ simulated events and is found to be 1.2 ± 0.2 compared to a fiat extrapolation 
(where the error is due to Monte Carlo statistics). 

Due to non-Gaussian tails in the distribution, the control region contains a 
mixture of Kl — > 37r° and vr^vr^ events. The number of vr^vr^ events is estimated and 
subtracted by using the distribution shape from tagged Ks events (after correcting for 
the contamination of accidentally mistagged Kl) and normalising this shape to the Kl 
distribution at small x^- 

Using this method, the background to the weighted Kl sample is found to be 
(8 ± 2) X 10~^ averaged over energy, where the error comes mostly from the Monte Carlo 
statistical error on the extrapolation factor. Changes in the choice of the control region lead 
to negligible variations in the estimated background. This procedure is done separately 
in each energy bin for the R computation, and the background is found to increase from 
^ at 70 GeV to ^ 0.2 % at 170 GeV. 




Figure 3: Distribution of for weighted Kl candidates (black triangles) compared to the 
shape for pure 27r° events derived from Ks candidates, normalized in the first bin. The 
excess of events in the Kl candidates in the region 36 < < 135 is used to extrapolate 
the background in the signal region (x^ < 13.5). 
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5.3 7r+7r~ background 

Semileptonic decays of Kl are selected by identifying the lepton as an electron (from 
E/p 1) or a muon (from associated hits in the muon counters and low E/p). For each 
kaon energy bin the event distribution in the {rriTrT^ , pt'^) plane is determined. The sample 
of Kl — >■ 7r"'"7r~ is then fitted as the sum of the contribution from true vr+vr" decays 
(determined from the shape of Ks tt+tt"), and from the Kg^ and K^^ distributions. 
Systematic uncertainties are evaluated determining the background in different control 
regions (contained in the domain 475 < m^jr < 505 MeV/c^ and p/^ < 2 x 10~^ GeV^/c^) 
and verifying the stability of the result for different sub-samples of data. The background 
in the signal region is (23 ± 2(stat) ± 4(syst)) x 10"^, where the result is averaged over 
the kaon energy spectrum. The larger contribution is from Ke3 decays. Figure ^ shows 
the event distribution in the variable p/"^, for Kl and Ks events, together with the result 
of the fit to the background. 

Kaon decays to 7r"'"7r~7 have been investigated and their effect on R is negligible. 

A small fraction (~ 0.5 x 10~^) of i^' — n~^7r~ are found in the Kl beam at large 
values of p/^, extending beyond (3-5) x 10"'^ GeV^/c^. These are due to scattering and 
regeneration on the collimators. Since these events are removed from the ir^ir^ sample by 
the pt'^ cut, but satisfy the centre of gravity cut and are kept in the 7r°7r° sample, they 
cause a correction to the observed value of the double ratio equal to (—12 ± 3) x 10~^. 
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Figure 4: Distribution of p/^ for Kl candidates fulfilling the 7r~^TT~ invariant mass cut and 
result of the fit for the background determination in the control region, including a signal 
component (derived from Ks events), a K^^ component and a i^^s component. The final 
signal region is p/^ < 2 x 10""^ GeV^. 



5.4 Energy and distance scale calibrations, stability, linearity 

The measurements of the kaon energy, the decay vertex and the proper time in 
the 7r°7r° mode rely on the measurement of the photon energies and positions in the 
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calorimeter. The calorimeter performances are studied and adjusted using the tt^tt^ events 
themselves, Kes decays where the electron energy measured in the calorimeter can be 
compared to the momentum measured in the spectrometer, and a so-called t] run which 
was a special run with a tt" beam striking a thin target near the AKS counter and 
producing vr" and r) (— >■ 77, Svr^) with known decay position. 

The overall energy and distance scale is adjusted using Ks ~^ tt^tt^ events. The 
beginning of the decay region is defined by using the AKS as veto. Therefore the distri- 
bution of the reconstructed vertex position exhibits a sharp rising edge at the nominal 
anti-counter position, as shown in figure The energy scale is set such that the aver- 
age value of the reconstructed decay position in a range centred around the anti-counter 
matches the value found in a Monte Carlo with the nominal scale. This procedure takes 
into account non-Gaussian tails in the energy resolution. The energy scale was checked 
daily and was stable within ±5 x 10~^ during the entire run. This measurement of the 
energy scale is checked using the data from the rj run. Another check is done using Tr^Tr*^ 
Dalitz decays and comparing the neutral vertex position with the vertex reconstructed 
using the charged tracks. The residual uncertainty on the overall energy scale is estimated 
to be ± 5 X 10~^ and this leads to an uncertainty of ± 5 x 10~^ on the double ratio. 

The procedure above assumes that the transverse size of the calorimeter is known. 
The size of the calorimeter is checked using K^s decays, comparing the reconstructed 
cluster position with the electron track impact point extrapolated to the calorimeter. 
This comparison is performed with an accuracy of 0.3 mm/m on the transverse scale, 
which leads to an uncertainty on the double ratio of ±3 x 10~^. 

The uniformity of the calorimeter response over it's surface is checked and opti- 
mised using Ke3 decays and tt" decays from the 77 runs. Residual non-uniformities in the 
calorimeter response over large distances are at the level of ~ 10"'^ and the corresponding 
systematic uncertainty on the double ratio is ±3 x 10~^. 

The calorimeter energy non-linearity is studied using Kg^ events, K^ s ~^ ti^tt^ 
events, and tt^ and rj decays from the 77 runs. Non linearities are found to be at the level 
of 0.3 % in the energy response between 3 and 100 GeV with possible energy offsets of 
± 20 MeV. These effects are taken as systematic uncertainties and lead to an uncertainty 
on the double ratio of ± 9 x lO"'^. 

Other uncertainties in the photon energy corrections (energy sharing between clus- 
ters and space charge effect), and from residual effects from the difference between the 
photon direction and the projectivity direction of the calorimeter are estimated to be 
± 4 X 10~^ on the double ratio. 

Adding all the previous effects together in quadrature, the total systematic uncer- 
tainty coming from the measurement of photon energies and positions is found to be 
± 12 X 10~^ on the double ratio. 

For charged decays, the vertex is measured from the reconstructed tracks and is 
completely determined by the detector geometry. As a check, the anticounter position is 
fitted on the Ks — ^ vr+vr" events (figure P^), and agrees to within 1 cm of the expected 
position. Uncertainties in the geometry and on the residual effect of the magnetic field in 
the decay region lead to an uncertainty on the double ratio of ±5 x 10~^. 

5.5 Accidental effects 

The effects of accidental activity in the detector are minimized by the simultaneous 
collection of data in the four channels. Most of the accidental activity is related to Ki 
decays and to muons originating from the production target. 
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Figure 5: Distribution of the reconstructed decay vertex in Ks events for (a) the vr^vr^ mode 
and (b) the tt+tt" mode. The rising edge corresponds to the position of the anticounter 
(represented by the dotted hne). 

The accidental rate can be measured directly from the activity in the detector within 
the readout window before each event. For example, the number of events affected by over- 
flows is the same for Kl — »■ tt+tt" and Kg tt+tt" to within the statistical error of 0.1 %. 
The correction to R due to this particular class of accidentals which prevent the event 
reconstruction in the drift chambers is negligible after removing events with overflows 
also from the samples of ir^ir^ decays, as discussed in section 4.4. Other accidentals may 
however cause a migration of events over the acceptance boundaries. This effect is also 
expected to be the same for Ks and Kl events due to their kinematical and geometrical 
similarity. 

A method to measure these effects is to overlay events with randomly triggered 
events which represent accidental activity in the beam, thereby artificially doubling the 
accidental intensity. The comparison of original and overlaid events provides a direct 
measurement of intensity effects in terms of gains and losses of events which are at the 
level of a few percent. As expected the effects tend to cancel in the double ratio. The 
correction due to accidentals is (—2 ± 14) x 10~^, where the error is dominated by the 
statistical uncertainty in the sample of the events used for the overlay analysis. 

The overlay technique does not cope with any additional detector activity in the 
Ks beam generated by the same proton which produced the Ks event. Studies of this 
background show that the effect on the double ratio can be neglected. 

5.6 Acceptance correction 

The Ks and Kl acceptances are made very similar in both modes by using the 
technique of weighting Kl events according to their proper decay time so that the effective 
longitudinal vertex distribution of Kl events is made equal to the Ks distribution. The 
weighting factor includes small contributions from interference and Ks decays in the Kl 
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beam, both of which are significant only for kaon energies above 140 GeV. The residual 
acceptance differences are studied with a full simulation of the beams and detector. The 
Monte Carlo sample is five times larger than the event sample. The main sources of 
systematic error in the acceptance correction are due to uncertainties in the positions 
and divergences of the beams and a significant component is caused by the dead column 
in the e.m. calorimeter. The overall tc^tt" reconstruction efficiency has been found to be 
the same for Ks and Kl decays using a Monte Carlo computation which includes the 
measured efficiencies of each drift chamber wire. The acceptance correction is evaluated 
and applied in each energy bin. The overall correction to the double ratio, averaged over 
kaon energies in the interval 70—170 GeV is (29 ± ll(MCstat) ± 6(syst)) x 10"^ 

6 Result 

The effects of the corrections applied to the double ratio are listed in table ^. The 
systematic errors shown in the first four lines of table § are dominated by the statistics 
of the control sample used in each study and are therefore statistical in nature. 

Table 2: Corrections applied to R and systematic uncertainties, in 10~^ units. 



Tagging Auls, ^^sl 


+ 18 


± 


11 


7r+7r^ trigger efficiency 


+ 9 


± 


23 


Acceptance 


+ 29 


± 


12 


Accidental effects 


- 2 


± 


14 


7r°7r° background 


- 8 


± 


2 


7r"^7r~ background 


+ 23 


± 


4 


Scattering/ regeneration 


- 12 


± 


3 


Energy scale and linearity 




± 


12 


Charged vertex 




± 


5 


Total correction 


+ 57 


± 


35 



The values of the double ratio in each bin of kaon energy are shown in figure |^, after 
trigger efficiency, tagging, background and acceptance corrections are included in each 
bin. The bins which are used in the analysis, between 70 and 170 GeV (decided prior to 
running the experiment), are shown with black symbols. The of the average is 25.7 (19 
d.o.f.). Extensive checks have been made to exclude systematic biases which could lead 
to a variation of R with the kaon energy as the data might suggest. As a further check, 
the double ratio R was computed in three additional energy bins, as indicated in figure 
These three extra points strongly disfavour the hypothesis of a linear trend in the data. 

The final result for the double ratio is i? = 0.9889 + 0.0027 + 0.0035, where the first 
error is statistical and the second is systematic. The corresponding value for the direct 
CP violating parameter is: 

Re {e'/e) = (18.5 + 4.5 + 5.8) x 10"^ 

Combining the two errors in quadrature, the result is: 

Re (e'/e) = (18.5 + 7.3) x 10"^ 

In conclusion, this new measurement of Re {e'/ e) confirms that direct CP violation occurs 
in neutral kaon decays. 
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Figure 6: Measured double ratio in energy bins. The points used for the measurement of 
Re {s'/e) are shown in black (energy range 70-170 GeV). 
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